Multiferroicity in rare-earth nickelates RNi0 3 
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We show that charge ordered rare-earth nickelates of the type RMO3 (R= Ho, Lu, Pr and Nd) are 
multiferroic with very large magnetically induced ferroelectric (FE) polarizations. This we determine 
from first principles electronic structure calculations. The emerging FE polarization is directly tied 
to the long-standing puzzle of which kind of magnetic ordering is present in this class of materials: 
its direction and size indicate the type of ground- state spin configuration that is realized. Vice 
versa, the small energy differences between the different magnetic orderings suggest that a chosen 
magnetic ordering can be stabilized by cooling the system in presence of an electric field. 

PACS numbers: 71.45. Gm, 71.10.Ca, 71.10.-w, 73.21.-b 



Introduction Complex oxides with simultaneous mag- 
netic and ferroelectric (FE) ordering -multiferroics- are 
attracting enormous scientific interest [TJ [2] . They offer 
the potential to control the magnetic order parameter 
by the FE one and vice versa -a very desirable prop- 
erty from a technological point of view [3]. Even if in 
quite a number of transition metal oxides both ferro- 
electricity and magnetism are present, magnetically in- 
duced FE polarizations observed so far are typically very 
small [4j |5j |6j [7]. This smallness is particularly pro- 
nounced in materials where multiferroicity relies on rela- 
tivistic spin-orbit coupling, which is intrinsically weak pQ . 

A few years ago it was pointed out that, theoretically 
at least, materials that are simultaneously magnetic and 
charge ordered can be multiferroic and potentially have 
a very large polarization [HJ^. To become multiferroic, 
however, an insulting oxide needs to meet an additional 
requirement: its symmetry has to be such that magnetic 
ordering can push a charge ordering pattern from site- 
centered towards bond-centered [8]. A large polarization 
results if the oxide is in addition electronically soft, so 
that inside it charge can easily be displaced. 

Here we show that precisely this scenario materializes 
in perovskite nickelates RNi03, where R is a rare earth 
element such as Ho, Lu, Pr or Nd. Consequently these 
nickelates can exhibit magnetically induced FE polariza- 
tions (P) that are very large, up to 10 \iC jew? . Such 
a polarization is two orders of magnitude larger than 
the one of typical multiferroics such as TbMn03[4] or 
TbMn205 [5 j. Also a very interesting fundamental point 
is associated with the symmetry of P in the rare-earth 
nickelates. To appreciate this aspect we have to bear 
in mind that in spite of their apparently simple chemical 
formula, the rare-earth nickelates are very complex mate- 
rials. They show an intriguing and only partially under- 
stood transition from a high temperature metallic phase 



into a low temperature insulating one. The nature of 
magnetic order in this low temperature insulating phase 
has been a long standing puzzle. Three different mag- 
netic structures have been proposed, two of which are 
collinear and one non-collinear, and so far experiments 
have not been able to differentiate between them. 

We show that all the proposed magnetic structures of 
the rare-earth nickelates are similar in the sense that all 
are multiferroic and very close in energy. However, differ- 
ent magnetic symmetries leave an individual fingerprint 
on the size and, in particular the direction of ferroelectric 
polarization P. In one type of collinear magnetic order- 
ing, for instance, P is parallel to the crystallographic 
6-axis; in the other it is perpendicular to this axis. With 
this theoretical result in hand, an experimental determi- 
nation of the direction and magnitude of P will reveal 
the precise type of magnetic ordering that is realized 
in the rare-earth nickelates. This fundamental observa- 
tion also suggests a practical application: in these nicke- 
lates it allows to control the realization of different mag- 
netic phases by cooling the material through its magnetic 
phase transition in an externally applied electric field. 

Lattice and Charge Order The metal-insulator tran- 
sition in RMO3 (R=Pr, Nd, Sm, Ho, Lu etc.) takes 
place at relatively high temperatures: Tmi =130 K (Pr), 
200 K (Nd), 400 K (Sm), 580 K (Ho) and 600 K (Lu) 
and is believed to coincide with the appearance of charge 
ordering and a simultaneous transition of the crystallo- 
graphic symmetry from orthorhomic Pbnm to monoclinic 
P2i/n [TOj [n]. The charge ordering is characterized 
by a nickel charge disproportionation of formally Ni 3+ 
into Ni( 3+5 ) + and Ni^ 3_5 ^ + , which form a simple two- 
sublattice, rocksalt-like, superstructure, see Fig. [I] To 
be specific we will consider in the following four repre- 
sentative members from the nickelate series: R = Ho, Lu, 
Pr and Nd. Of these R=Ho and Lu are small rare-earth 
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FIG. 1: (Color online) Schematic view of the charge and mag- 
netic ordered a) S-type b) T-type and c) N-type non-collinear 
magnetic structure of RNiOs. Ni 2+ and Ni 4+ are shown in 
red and orange. Arrows represent Ni spins. In a) and b) the 
zig-zag spin chains are indicated. In c) the FM planes per- 
pendicular to [101] are highlighted. The vector P shows the 
direction of polarization obtained in our calculations. 



ions, with a magnetic ordering temperature Tn below the 
metal-insulator transition: Tn < Tmi (for Ho Tn=145 K, 
for Lu 130 K). The nickelates with the larger rare-earth 
ions, Pr or Nd, are different in that their metal- insulator 
transition coincides with the appearance of magnetic or- 
dering: Tn =Tmi- The monoclinic P2i/n crystal struc- 
ture contains two inequivalent Ni positions and three 
inequivalent oxygen atoms (Oi, O2, O3) [TQl [TTJ [T2] . 
Experimentally the charge ordering reflects itself in an 
oxygen breathing distortion of the Ni06 octahedra and 
induces different magnetic moments on the two inequiv- 
alent Ni atoms (1.4/1.4 \±b and 0.6/0.7 \±b for Ho/Lu, 
respectively) [T3l H4]. 

Magnetic Order The magnetic neutron diffraction 
peaks are characterized by the propagation vector Q = 
( ^ , 0, \) p~2j [IHl [15] . In this class of nickelates the exper- 
iments suggested three possible magnetic structures. In 
the following we label two collinear magnetic structures 
as S, T , and the non-collinear one as N (see Fig. [I]). In 
both the S and T-type magnetic structure the Ni spins 
within the afr-plane form ferromagnetic zig-zag chains. 
Adjacent zig-zag chains in the same plane are coupled 
antiferromagnetically. The stacking of the zig-zag chains 
along the c-axis differentiates between S and T-type or- 
dering: in S-type (Fig. |TJl) the zig-zag spin-chains are 
all pointing in the same direction, whereas for T-type 
(Fig. [TJd) the zig-zag chains in adjacent planes are point- 
ing in alternate -fa and -a directions. This implies for 
T-type that the spins in planes perpendicular to the [- 
1,0,1] direction are pointing in the same direction, see 
Fig. [TJd. Along this [-1,0,1] direction these ferromagnetic 
planes are ordered in a TTIITTII fashion [9]. In the non- 



collinear N-structure (Fig. [TJ: [16 ) all spins lie in the ac 
plane, with spins in a plane perpendicular to [1,0,1] point- 
ing in the same direction. When moving from plane to 
plane along [1,0,1] the spins rotate within the ac-plane. 
The N-type magnetic structure therefore corresponds to 
a spin-spiral. 

Ab Initio Results For the crystal and magnetic struc- 
tures outlined above we performed a set of density func- 
tional calculations using the projector augmented- wave 
(PAW) method and a plane-wave basis sets as imple- 
mented in VASP [17 . We include the strong Coulomb 
interactions between the Ni 3d electrons, in SGGA+U 
[IHl [19J [20] calculations for U =8 eV and a Hund's rule ex- 
change of Jh = 0.88 eV. Starting from the experimental 
centrosymmetric crystal structures we compute the elec- 
tronic structure for T, S and N magnetic order. In all the 
calculations, the difference in total energy between the T 
and S type ordering is very small, within the numeri- 
cal accuracy. We evaluate the electronic contributions to 
the polarization with the Berry phase method within the 
PAW formalism [ST] [22]. In Table [j] we report the elec- 
tronic contributions to the FE polarization for the T and 
S type magnetic ordering. For the non-collinear N-type 
ordering these calculations are numerically extremely de- 
manding. We have therefore computed the polarization 
for H0MO3 in its experimental crystal structure only, re- 
sulting in P c =-110 nC/cm 2 and P a =20 nC/cm 2 . It is 
remarkable that for T-type the polarization is large, in 
the ac-plane and predominantly along the a-axis, for S- 
type it is large and strictly along the 6-axis and for N-type 
it is weak, in the ac-plane and mostly along c. 

S- and T-type Magnetism: Origin of Multiferroicity 
Having established that the magnetic ordering induces a 
very significant FE polarization in the rare-earth nicke- 
lates, we now clarify its microscopic origin and explain 
its different direction for T and S-type order. In the 
P2i/n structure the corner-sharing NiC>6 octahedra are 
distorted and tilted due to the so-called GdFeOa distor- 
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TABLE I: FE polarization P of RNi0 3 (R=Ho, Lu, Pr) in 
\iC I cm 2 for both the experimental centrosymmetric (Ref . [10] 
for Ho, Lu and Ref. [TI] for Pr) and the relaxed crystal 
structure, for either S or T-type magnetic ordering, with 
P2i and Pn symmetry, respectively. Atomic positions for 
NdNiOs taken from experiments in the Pbnm orthorhombic 
setting [23] and then relaxed without changing the Bravais 
lattice. 
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FIG. 2: (Color online) Schematic arrangement of oxygen- 
induced dipoles and resulting FE polarization for the P2i/n 
crystal structure, in the top layer of Fig. [TJl and b. 



tion, which causes a counter-rotation of neighboring oc- 
tahedra in the afr-plane, see Fig. [2] Neighboring planes 
along the c-axis display a similar in-plane oxygen dis- 
placement pattern. The end result of the distortions is 
that each Ni ion is surrounded by two crystallographi- 
cally inequivalent oxygen sites in the afr-plane: O2 and 
O3, see Fig. [2] Along one direction chains of 02-Ni-02-Ni 
bonds form and perpendicular to it, all bonds are 03-Ni- 
O3-M. In spite of these distortions the crystal structure 
is still centrosymmetric and the material therefore para- 
electric. However, the formation of zig-zag spin chains in 
the afr-plane breaks this inversion symmetry. 

The magnetic ordering along the M-O2 chain direc- 
tion can be denoted as Ni|-02-Ni|-02-Ni|. In this struc- 
ture the O2 sites become inequivalent as one O2 is in 
between charge disproportionate Ni ions with parallel 
spin and the other one between antiparallel Ni spins. 
Thus the oxygen sites split into O2 and O2' , see Fig. |2j 
Besides the charge disproportionation on the Ni sites, 
now also the inequivalent oxygen atoms O2 and 02' 
charge-polarize. For the Ni-03 chain the situation is 
similar and the splitting is into O3 and O3/. The in- 
equivalence of the four in-plane oxygens ions surround- 
ing a Ni is directly reflected by their different Born ef- 
fective charges Z*: when for instance R=Lu we find 
Z*(0 2 >) = -1.25e, Z*(0 2 ) = -6.37e, Z*(0 3 >) = -1.59e 
and Z*(Os) = — 1.73e. For other rare-earths we observe 
similar trends. The resulting inequivalence of oxygen ions 
situated on the nickel bonds causes a partial shift away 
from a nickel site-centered charge ordering to a Ni-Ni 
bond centered charge ordering. The resulting net dipole 
moment of each nickel-oxide afr-plane has a finite projec- 
tion along both the a and 6-axis, see Fig. [2] This dipole 
formation is reminiscent of the mechanism for ferroelec- 
tricity proposed in Ref. [24], where it was discussed in 
the context of HoMn03. For S-type magnetic order the 
dipole moments of different nickel-oxide afr-planes add 
up to a net polarization along the 6-axis. The following 
symmetry argument underlies the fact the polarization 




FIG. 3: (Color online) Ferroelectric polarization after an el- 
ementary translation along the c-axis, corresponding to the 
second layer Fig. [I] for a) S- and b) T-type magnetic order- 
ing. 



vanishes in all other directions. Translating an ab plane 
along [001] interchanges the O2 and O3 positions, but 
leaves alternation pattern of parallel/anti-parallel spin 
bonds invariant, see Fig. |3^i. The O2-O3 interchange is 
equivalent to a rotation around the b- axis by tt. As in 
this case also P rotates by tt around 6-axis, the net po- 
larization points fully along 6, see Fig. [3^. 

By virtue of the same argument for T-type magnetic 
ordering the polarization along the b- axis vanishes, see 
Fig. [3)3. In this case a translation along [001] again inter- 
changes the O2 and O3 positions, but the bonds of paral- 
lel spins transform into antiparallel ones and vice versa. 
This implies that compared to the S-type structure, an 
elementary translation along [001] in T-type generates an 
additional reflection of P in the a-axis, see Fig. (3}d. Con- 
sequently the net in-plane polarization of the system is 
now along a. Via a similar argument one finds in T-type 
a finite polarization along the c-axis, which is forbidden 
in the S-type structure. 

Non-collinear phase: Origin of Multiferroicity The N- 
type structure corresponds to a magnetic spiral with 
propagation vector Q= [1/2,0,1/2] and rotation axis e = 
[0,1,0]: it is a spiral with spins rotating in the ac-plane. 
The spiral structure implies that in the direction of Q 
the angle between spins in successive planes is invariant, 
so that magnetic bonds are equivalent. The N-type mag- 
netic ordering can therefore not give rise to multiferroic- 
ity via a modulation of the charge ordering, as is the case 
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in the S- and T-phase. However, in spin-spirals the rela- 
tivistic spin-orbit coupling directly causes a FE polariza- 
tion pQ. Symmetry dictates that this polarization arises 
in the direction e x Q, corresponding to [1,0,-1] in our 
case. Indeed, our ab-initio calculations including spin- 
orbit coupling give P = (20,0,-110) nC/cm 2 , in full 
agreement with these symmetry considerations. The size 
of the polarization is comparable to that of TbMnOs [4] 
and thus much weaker than that for the S- and T-type 
ordering, which is a generic feature of systems where mul- 
tiferroicity is caused by spin-orbit coupling. 

Relaxed structure We have also computed the ionic con- 
tribution to P that arises from lattice distortions induced 
in the FE phase. Starting from the experimental cen- 
trosymmetric P2i/n crystal structure with a magnetic 
supercell we relax unit cell ionic positions, allowing for 
a lower symmetry structure to develop. Details of the 
structural relaxations will be published elsewhere [25] . 
The resulting total FE polarizations are reported in Ta- 
ble [l| For the S-type magnetic state the induced lattice 
distortions enhance P for all nickelates that we have con- 
sidered. For T-type also a reductions of the polarization 
occurs in some materials, but still P stays very large: 
- 10/iC/cra 2 for Ho/Lu, - 8/iC/cm 2 for Pr/Nd. For 
T-type the polarization along the a-axis, P a , dominates 
over P c . We find that P c is so small due to a partial 
cancelation of ionic and electronic contributions to P c , 
an effect that was also observed in HoM^Os 6j and in 
TbMn0 3 [26 . For the S-type structure we find total 
P b ~ 7/iC/cm 2 for Ho/Lu, and - 3/iC/cm 2 for Pr/Nd. 
The generic observation is that nickelates with small rare- 
earth ions (Ho and Lu) tend to show the largest polar- 
izations. 

Conclusions On the basis of theoretical calculations we 
predict that the perovskite nickelates RMO3 (R - rare 
earth) are multiferroic in their low-temperature insulat- 
ing magnetic phase. We show there are different mech- 
anisms at play for magnetically-induced ferroelectricity 
in RMO3. The S and T collinear magnetic spin con- 
figurations give a remarkably large polarization (~ 10 
/iC/cm 2 ) along the b axis and in the ac plane, respec- 
tively. It is driven by non-centrosymmetric arrangements 
of spins which force the charge-ordering to shift from site- 
centered to partially bond-centered. The estimated po- 
larization in the S and T collinear cases is much larger 
than that arising from a relativistic spin-orbit related 
mechanism, which is at play in the N-type spin-spiral 
state. An experimental determination of the direction 
and magnitude of P can therefore solve the long-standing 
puzzle related to the magnetic ground-state in nickelates. 
The fact that according to our calculations the energies 
of the S- and T-type magnetic structures are very close 
but the directions of polarization are quite different, sug- 
gests one can stabilize one or the other by cooling in an 
appropriate electric field, so that in effect an applied ex- 
ternal electric-field can control the realization of different 



magnetic phases in these nickelates. 
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